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Abstract 

Jets observed in hadron-hadron scattering contain a contribution from the 
"underlying event" that is produced by spectator interactions taking place in- 
coherently with the major parton-parton collision, due to the extended com- 
posite structure of the colliding hadrons. Using a recent measurement of the 
double parton interaction rate, we calculate that the underlying event may 
be 2-3 times stronger than generally assumed, as a result of semi-hard per- 
turbative multiple-parton interactions. This can have an important influence 
on the inclusive jet cross section at moderate values of Ex, persisting at the 
5 — 10% level to the largest observable Et- We show how the underlying 
event can be measured accurately using a generalization of the method first 
proposed by Marchesini and Webber. 
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I. INTRODUCTION 



The inclusive jet cross section da/dEx, averaged over a small range of pseudorapidity, 
is an important object for study because it tests perturbative QCD at the highest Q 2 scale 
currently possible Beyond its potential for detecting physics beyond the standard 

model or confirming QCD, the jet cross section is even beginning to play a role in the global 
data fitting used to measure parton distribution functions . It is therefore important to 
carefully consider all systematic effects that influence the interpretation of the measurement. 

Among those effects is the "underlying event" generated by spectator interactions that 
can occur concurrently with a major parton-parton collision, due to the extended composite 
structure of the colliding hadrons. Simulations such as the HERWIG Monte Carlo ]5| include 
a "soft" underlying event that is modeled by a parametrization of minimum bias data.[| 
However, there may also be a hard underlying event: particles can be created by incoherent 
parton-parton interactions at momentum transfers that are sufficiently large for perturbative 
QCD to be a useful approximation, but much smaller than that of the interaction mainly 
responsible for a given high Et jet event. 

Attempts have been made in the past to predict underlying event contributions, including 
the perturbative part that is the focus of our attention |7|-|9| . We make a new estimate here, 
based on a recent direct measurement [|l(J of the rate for double parton interactions in 
hard scattering. In the following sections we derive the prediction, show how it can be 
independently tested by experiment, and discuss its effect on the inclusive cross section. 



II. HARD UNDERLYING EVENT 

To generate typical semi-hard perturbative final states, HERWIG 5.8 [H] was used to sim- 
ulate QCD 2^2 hard scattering at the Tevatron energy yfs = 1.8 TeV, with the minimum 
p± parameter set to a fairly small value PTMIN = 2.0 GeV that is nevertheless large enough 
for comfort with the perturbative calculation. This choice yields a cross section a = 50.6 mb, 
which is approximately equal to the full inelastic non-diffractive cross section [|ll|] , so we ap- 
pear to be taking an extreme point of view in which that cross section is mainly generated 
by perturbative ("minijet") interactions. The point of view is actually not so extreme, since 
including s-channel unitarity effects, e.g., by an eikonal model, would substantially reduce 
the cross section on the basis that once one interaction has taken place, additional inter- 
actions do not really add to the inelastic cross section |Ti2| , |TE|] . The soft underlying event 



feature of HERWIG was turned off (PRSOF= 0), since our goal is to study the perturbative 
part of the underlying event. 

If o"i and <72 are cross sections for distinguishable rare parton-parton interactions, the 
cross section for both interactions to occur in the same event can be written as 

_ °" lCr 2 m 

^double — , UJ 

Ocff 



Minimum bias events provide at best an imperfect model of the underlying event, since for 
example they always contain particles; while the soft underlying event can sometimes be absent, 
as shown by the finite survival probability for inelastic events with large rapidity gaps || . 
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with the parameter a c s conveying the probability for double parton interaction. a e g carries 
nonperturbative information beyond the scope of ordinary parton distribution functions, 
since it relates to correlations between partons within a single hadron. It has recently been 
measured to be <7 e fj = 14.5 ± l.Tt^mb [jTOf . 

For a fixed small value of a±, corresponding to some rare hard scattering of type 1, 
Eq. (d) gives a-2/oeff as the probability for a second rare scattering of type 2. If this type 
2 scattering is not rare, and it occurs randomly, the obvious generalization of Eq. ([!]) is a 
Poisson distribution in the number n of type 2 interactions that accompany a given type 1 
interaction: 

P n = - x e- a (2) 
n! 

with mean number a = <7 2 /cr e f[. 

Our model for the hard underlying event therefore consists of superimposing n of the 
above minijet events from HERWIG, with a Poisson probability distribution (|2]) of mean 
a = 50.6 mb/14.5 mb = 3.49 . In doing this, we ignore energy-momentum conservation effects 
in the sense that the parton momentum distributions are taken to be the same as in single- 
interaction events. This assumption should be adequate because the minijet interactions 
don't require a great deal of energy, and thus come from partons at rather small x. (More 
elaborate models based on PYTHIA |TJ]], and || based on HERWIG, have incorporated 
these energy-momentum constraints.) Meanwhile, the Poisson assumption could actually 
underestimate the frequency of 3 or more interactions, since these may be enhanced by "hot 
spots" — e.g., by strong spatial correlations between the partons in a beam or target hadron 



associated with constituent quarks [15 



In this paper, we focus specifically the case of the underlying event in dijet production. 
However, it will be valuable to look also at W, Z°, or lepton pair production, where the 
hard scattering makes a color-neutral object, so there is no radiation associated with final 
state jets. These processes also share the advantage that the hard scattering is initiated by 
qq, which generates less radiation than qg or gg. Single jet final states such as W + 1 jet 
and 7 + 1 jet should also be studied. 

We generate events with jets of transverse energy Et > 100 GeV using HERWIG, and look 



at the total transverse energy inside a control region cone R = J (rj — r/ ) 2 + (<p — <f) ) 2 < 0.7, 
which is the size typically used to define jets in the pseudorapidity-azimuthal angle "Lego" 
plane. The underlying event is by definition uncorrelated with the jets, so to avoid the 
majority of the jet Et, the control cone is centered a distance R = 2.0 away from both jet 
axes, in a manner to be described fully in Sect. [TQ] . 

Fig. 1 shows the distribution of transverse energy Et in the control region cone from the 
HERWIG soft underlying event (dashed curve) and from the hard underlying event model 
(dotted curve). The average values, which are relevant for the correction to the inclusive jet 
cross section as will be discussed in Sect. [TV], are (Et) = 0.84 GeV for the soft background 
and (Et) = 1.90 GeV for the hard background. For comparison, Fig. 1 also shows the 
predicted contribution from radiation associated with the hard scattering (solid curve), as 
given by HERWIG with the underlying event turned off. It has (Et) = 2.37 GeV. 

The hard underlying event model predicts a substantially larger (Et) than the HERWIG 
soft model. Taken as estimates of the average underlying event Et that should be added 
to a perturbative QCD calculation to predict the observed inclusive jet cross section, the 
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difference is significant as will be discussed in Sect. [TV] - the more so if both mechanisms 
operate concurrently, so their contributions to E T should be added. 

Before turning to the influence on inclusive jet cross sections, we first consider how to 
measure the background event better. 



III. OBSERVING THE UNDERLYING EVENT 

An obvious way to measure the background event contribution to jets is to look at the 
E T distribution in a "control region" cone of the same radius that is used to identify jets, 
as discussed in the preceeding Section. Fig. 2 shows the predicted probability distributions, 
which correspond to adding the contributions described by Fig. 1: HERWIG with its soft 
underlying event model included (dashed curve), HERWIG with its soft underlying event 
portion replaced by the hard underlying event model of Sect. || (dotted curve), and HERWIG 
with both underlying event models operative (solid curve). All three curves share contribu- 
tions from the substantial QCD radiation generated by the principal hard scattering and its 
color connections to the beam particles. This makes the curves somewhat similar, but the 
differences are large enough that a measurement of this control Et distribution would give 
a useful indication of the background event level. 

To measure the underlying event more accurately, we generalize a technique first advo- 



cated by Marchesini and Webber |16| - which has apparently not yet been tried experi- 
mentally in its original form. The essence of the technique is to study Et simultaneously 
in two regions that are separate from each other and separate from any jets that define the 
final state under study. 

We describe the proposed generalization in the context of our application to dijet events. 
Let (771, 0i) and (772,02) be the locations of the jets in the Lego plane. For this study, we 
make the approximation that the jets are back-to-back in azimuth: |0i — 02 1 = 7r. Label the 
jets so that r/i < 7/ 2 and let 



K 



y/4/[n 2 + (Vi-m) 2 ] -1/4, (3) 



rejecting the small fraction (5.6%) of events that have \r]i — 1] 2 \ > 2.48, for which k is 
undefined. Then consider the two points (r] a ,(j) a ) and (r]b,4>b) where 

Va = Vb = {Vi + ?72)/2 + kit (4) 
= 0i ± (tt/2 - k\t] X - 772!) . (5) 

These points are at the same pseudorapidity r\ a = r] b , and are well separated from each other 
in azimuth by |0 a — 0{,| = it — 2K,\r] 1 — 772 1 which varies from 2.28 to 3.14 . Both points are 
separated from both of the jet axes by a distance of exactly 2.0 in the Lego plane. (The jet 
that lies farther from them in 77 is therefore closer to them in 0.) We can define a similar pair 
of points by letting n — > —k in Eq. (|j) and 0i — > 2 in Eq. (|[). For the present purpose, we 
use only the pair of points with r\ a = i] b closer to 77 = to reduce any kinematic suppression. 

Let S a and S b denote the total transverse energy E T deposited in cones of radius 0.7 
centered on these two points. Fig. 3 shows some typical locations of the cones with respect 
to the jet axes. From the HERWIG simulation (taking the jet axes from the directions of their 
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partons, without modification by initial state radiation for simplicity), we find ± 0.64 for 
the mean and standard deviation in rj and 2.59 ± 0.26 for the mean and standard deviation 
in \4> a — (pb\- (In practice, it may be better to keep the control cones centered at fixed values 
of rj, at the expense of letting their distances from the jets vary somewhat, to avoid effects 
due to the ^-dependence of detector corrections.) 

The cones a and b are those already used in Sect. |IJ to study the underlying event 
background. The two cones are equivalent, and only one was used for each event. Fig. 2 
can therefore be interpreted as the probability distribution for S a or equivalently for S&. 

Because there is no intrinsic difference between the two control cones, averaging over 
events would give (S a ) = (Sb). It follows that 

(Sa) = (S b ) = (^y^) = (Sm) + (Smin) (6) 

where 

Sdifr — \S a — Sb\/2 (7) 

S min = min(S a , S b ) . (8) 



This separation into S'difr and S m i n is useful, as pointed out in Ref. jig] , because next-to- 
leading order perturbative corrections to the principal hard scattering contribute in at most 
one of the two regions, and hence contribute only to S'difr; while the underlying event (like 



minimum bias events |17]]) is expected to have positive correlations over long distances in 



(rj, </>), so its contribution to Sdifr is suppressed while its contribution to S m i n is strong. Thus 
measuring the distributions of Sdifr and S m i n separately will be much more revealing than 
the distribution of S a or S b of Fig. 2 alone. 

Fig. 4 shows the predicted probability distributions for Sdifr- As anticipated above, Sdifr 
is not very sensitive to the choice of model for the underlying event. Testing this distribution 
against experiment would therefore be a good way to test QCD in a manner that is not very 
sensitive to underlying event effects. 

Fig. 5 shows predicted probability distributions for S min . As anticipated above, S min is 
very sensitive to the underlying event model. Testing this distribution against experiment 
would therefore be an excellent way to measure the strength of the underlying event. 

The method of Marchesini and Webber [K| could be generalized further in the hunt for 



the background. For example, the transverse energy could be measured in all four of the 
control regions that can be defined by R = 0.7 cones centered 2.0 units in (rj, 0) from both 
jets. These regions are shown in Fig. 3 for jets that are back-to-back in 0. The distribution 
of the minimum of the four transverse energies, or the distribution of the sum of the two 
smaller ones, would be especially sensitive to the underlying event. Other cone sizes would 
also be of interest, but R = 0.7 applies most directly to the correction to the inclusive jet 
cross section. 



IV. CORRECTION TO PQCD INCLUSIVE JET CROSS SECTION 

The influence of the underlying event on the inclusive jet cross section is calculated as 
follows. Let 
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F(Et) EE ^ (9) 

be the single hard scattering contribution to the jet cross section, which is predicted by 
standard pQCD techniques on the basis of parton distribution functions ||. Let 

G (Et) . §r (10) 

be the probability distribution for additional E T inside the jet cone contributed by the 
underlying event, normalized to G(E) dE = 1. The observable jet Et is the total of the 
two contributions, so the observable inclusive cross section is given by 

F ohs (E T ) = ^ = f°° dE, G(E 1 ) [°° dE 2 F{E 2 ) 5{E X + E 2 - E T ) 
a E T Jo Jo 

= [ ET G(Ei)F(Et- E 1 )dE 1 . (11) 
Jo 



It is assumed here that the underlying event contribution is small enough not to significantly 
shift the apparent jet axis. In practice there would be a small additional increase in the 
average jet energy due to the jet-finding algorithm's tendency to maximize the Et included 
in the jet. 

Fig. 6 shows the fractional increase in the inclusive jet cross section caused by the HERWIG 
soft background event (dashed curve); by the hard underlying event model of Sect. [TT] (dotted 
curve); or by including both underlying event models (dot-dash curve). The effect of the 
underlying event is rather large at the lower values of jet Ex, so it must be measured rather 
well before a meaningful comparison can be made between the observed jet cross section and 
the pQCD prediction. The effect of the underlying event is significant even at the highest 
jet Et shown, where the cross section is raised by approximately 3%, 6%, or 9% under these 
three assumptions. The fractional increase caused by adding the hard underlying event to 
the prediction with the soft underlying event already included is shown by the solid curve in 
Fig. 6. It is nearly identical to the the increase caused by adding it to the pQCD prediction 
alone (dotted curve). 

The procedure used by experimenters |],|] to take account of the underlying event, along 
with background from simultaneous events ("pile- up"), is to subtract a constant value from 
each measured jet Et- We examine the accuracy of that approach next. 

An average underlying event contribution U can be defined by the exact equation 

F ohs {E T ) = F{E T -U) . (12) 

Since G(E) falls rapidly with E, it is natural to make a linear expansion of F(E) in the 
neighborhood of Et in Eq. (0). The underlying event contribution is then given by 

POO 

U^E= G(E) E dE , (13) 
Jo 

which corresponds to the experimental procedure of approximating the background contri- 
bution in each event by the average value. This average is subtracted from each measured jet 
Et before obtaining the inclusive cross section that is compared with pQCD. The value used 
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is 0.9 — 1.1 GeV, which is similar to the HERWIG soft background prediction (Ex) = 0.84 GeV 
found above. It corresponds to an underlying event level of ~ 0.6 — 0.7 GeV per unit area 
in (r/,0). 

I find that the linear approximation works rather well for the models considered here, even 
in the most extreme case where both soft and hard background are included. For example, 
in that case (Ex) = 0.84 + 1.90 = 2.74 GeV, while the true effect corresponds to U — 3.0 for 
jets of Et = 100 GeV. However, in real experiments, there is a further background from true 
minimum bias collisions that occur between other pp pairs at high luminosity. This can easily 
raise the background level to the point where the linear approximation breaks down. In that 
case, the CDF analysis method |I[ of parametrizing F(Ex) and fitting the parameters to the 
experimental results, which is useful to unfold other detector effects anyway, can include this 
one as well. The essential need is to allow for the event-to-event fluctuations in background. 

In the linear approximation, the fractional effect on the inclusive cross section can be 
written in the form 

F ohs (E T ) - F(E T ) ^ nU_ 
F(Et) Ex 

where n is the local effective power law defined by F oc E^ n , i.e. 

n(E T ) = -d(]nF)/d(\nE T ) . (15) 

Over the range 50 GeV < E T < 400 GeV, n(E T ) rises from ~ 5.5 to ~ 12.5. Its large 
value, which represents the rapid fall of the inclusive cross section with Ex, enhances the 



effect of the background according to Eq. ([14]) , as has been emphasized recently [ 18| . For 
example, it implies that an additional 1-2 GeV of background Ex, which may be present by 
the mechanism of Sect. ||, will raise the inclusive cross section by 3-6% at Ex = 400 GeV, 
which is consistent with Fig. 6. 



V. CONCLUSIONS 

The point of this paper is that the inclusive jet cross section may contain an underlying 
event contribution that is 1-2 GeV larger than generally assumed, due to incoherent multiple 
semi-hard interactions that accompany the hard scattering. The calculation presented here 
is less elaborate than previous methods of assessing these multiparton interactions |7|-|9| . 
but it has the advantage of being constrained by a recent measurement of the double parton 
interaction rate 11(11 . 



In view of the importance of the inclusive jet cross section, it is urgent to settle the 
question of underlying event level definitively by measuring it as described in Sect. 
This will also serve as a test of the assumptions used to make the estimate, which involve 
interesting unexplored areas of non-perturbative QCD. 

We have focused on Ex in an R = 0.7 cone, because that is the relevant quantity 
for estimating the incoherent background under a jet. The background will also have an 
important influence on "jet shape" observables such as the cone size dependence ip(r) = 
J2n<r E T / J2n<R E T 0- To study the underlying event mechanism in more detail, it would 
be useful to apply the technique based on Marchesini and Webber [fTH] that is discussed in 
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Sect, m using regions of different sizes as well. It would also be useful to apply the technique 
to different hard processes such as W and W + ljet production, and for comparison to 
minimum bias events. 
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FIG. 1. Normalized probability distributions for Et in R = 0.7 control cone. Solid curve 
= HERWIG with soft underlying event off; dashed curve = soft underlying event from HERWIG; 
dotted curve = hard underlying event from Sect. [fij. 
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FIG. 2. Predicted probability distributions for Et in R = 0.7 control cone. Dashed curve 
= HERWIG including soft underlying event; dotted curve = HERWIG with soft underlying event 
replaced by hard underlying event model of Sect. O; solid curve = HERWIG with both soft and 
hard underlying event included. 
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FIG. 3. Typical "control" regions a,b,a',b' for measuring the transverse energies S a and Sb- 
Configurations shown are for jets J\ and J2 with \rji — 772 1 = 0.0, 1.0, 2.0. The centers of the control 
regions are a distance 2.0 from both jet axes. The distance from a to b is the same as the distance 
from a' to b', although this is not apparent because the azimuthal cylinder is cut for display at 
\<fi\ = 7T. We use the pair ab or a'b' (here ab) that is closer to i] = 0. 
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FIG. 4. Normalized probability distributions for Sdis = \S a — <S&|/2, where S a and S& are 
the total transverse energies in the two R = 0.7 control regions. Solid curve = HERWIG with 
soft underlying event off; dashed curve = HERWIG with soft underlying event on; dotted curve 
= HERWIG with soft underlying event replaced by hard underlying event from Sect. ||; dot-dash 
curve = HERWIG with both soft and hard underlying event. 
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FIG. 5. Normalized probability distributions for S m i n = min(S' a , Sb)- As in Fig. 4, Solid curve 
= HERWIG with soft underlying event off; dashed curve = HERWIG with soft underlying event 
on; dotted curve = HERWIG with soft underlying event replaced by hard underlying event from 
Sect. ||; dot-dash curve = HERWIG with both soft and hard underlying event. 
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FIG. 6. Fractional increase in the inclusive jet cross section produced by adding HERWIG 
soft underlying event (dashed curve); the hard underlying event model of Sect. [0] (dotted curve); 
or both underlying event contributions (dot-dash curve) to the pQCD prediction. The fractional 
increase produced by adding the hard underlying event to the soft+pQCD prediction is shown by 
the solid curve. 
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